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Abstract 

We construct a linear sigma model for open-strings ending on special 
Lagrangian cycles of a Calabi-Yau manifold. We illustrate the construction 
for the cases considered by Aganagic and Vafa(AV). This leads naturally 
to concrete models for the moduli space of open-string instantons. These 
instanton moduli spaces can be seen to be intimately related to certain 
auxiliary boundary toric varieties. By considering the relevant Gelfand- 
Kapranov-Zelevinsky (GKZ) differential equations of the boundary toric 
variety, we obtain the contributions to the worldvolume superpotential 
on the A-branes from open-string instantons. By using an ansatz due to 
Aganagic, Klemm and Vafa (AKV), we obtain the relevant change of vari- 
ables from the linear sigma model to the non-linear sigma model variables 
- the open-string mirror map. Using this mirror map, we obtain results 
in agreement with those of AV and AKV for the counting of holomorphic 
disc instantons. 
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1 Introduction 



D-branes have come to play a fundamental role in our understanding of string 
theory, especially for the case Calabi-Yau (CY) compactifications. Subsequent 
to the initial description of D-branes wrapping supersymmetric cycles of CY 
manifolds Jl|, it was realised that A and B-type D-branes^ are closely related 
to the objects that are relevant to Kontsevich's homological mirror symmetry 
proposal The fact that mirror symmetry, generally speaking, exchanges 

A-branes and B-branes has begun to play a crucial role in understanding several 
aspects of stringy geometry on Calabi-Yau manifolds. 

One of the applications of mirror symmetry in type II compactifications has 
been the computation of closed string instanton corrections to various physical 
quantities. With the inclusion of D-brane sectors for strings moving on Calabi- 
Yau backgrounds, one has, in addition, to consider a generalisation of mirror 
symmetry that also includes the open string sectors. In particular, one may 
expect that mirror symmetry in the open and closed string sectors taken together 
is relevant to understanding the contribution of open-string instanton corrections. 

In contrast to the case of purely closed strings on CY, it has been pointed 
out that various quantities in the world-volume theories of both A and B type 
D-branes acquire open-string instanton corrections. In general it is expected that 
the F-terms in the world-volume theory have open-string instanton corrections for 
A-branes whereas the D-terms are the ones that acquire instanton corrections for 
B-branes Various general considerations relevant to open-string instanton 

corrections in the case of compact Calabi-Yau manifolds as well as non-compact 
Calabi-Yau manifolds have been studied earlier ||[7|. 

However the explicit calculation of open-string instanton corrections to the 
world-volume superpotential for A-branes using mirror symmetry have first ap- 
peared recently in the important work of Aganagic and Vafa || and Aganagic, 
Klemm and Vafa ||. The authors in these papers have computed the world- 
volume superpotential for some non-compact A-branes in non-compact Calabi- 
Yau manifolds (in particular, those that can be described without a world-sheet 
superpotential), by appealing to mirror symmetry and using the holomorphic 

The terminology arising from the topological theory that support these branes, in accor- 
dance with the open-string world sheet boundary conditions that preserve some amount of 
supersymmetry in either case. We shall refer to the D-branes as A-branes/B-branes indicating 
the appropriate boundary condition. 
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Chern-Simons action on the B-brane side to reduce the problem to a more 
tractable computation. This computation uses and extends the techniques first 
described in [ID, IT|, whereby the mirror transformation, in the presence of D- 
branes, can be explicitly implemented as a world-sheet duality transformation 
involving the fields of a Gauged Linear Sigma Model (GLSM) that describes 
these D-brane configurations. Remarkably, the authors of |8|,|9| showed that the 
(instanton generated) A-brane superpotential thus computed conformed exactly 
to the strong integrality predictions of such quantities from completely different 
considerations of large N duality |~2] . 

We briefly recall here the manner in which this integrality predictions were 
first made. In [13] it was conjectured that the large N limit of Chern-Simons the- 
ory on S 3 is dual to A-type closed topological string theory on the resolved coni- 
fold, with a particular identification between the parameters of the two theories. 
This conjecture was verified at the level of the partition function on both sides. 
The conjecture was extended to the observables of the Chern-Simons theory, 
namely the knot invariants in [[HJ] . By identifying the extra D-branes (associated 
with the knots on the Chern-Simons side) on the resolved conifold, non-trivial 
predictions were made in [HJ] regarding the structure of the topological string 
amplitudes, using ideas of |0|. (Knot theoretic aspects of the conjecture were 
verified in [[14]].) Further, from the fact that topological disc amplitudes compute 
superpotential terms in iV = 1 d = 4 gauge theories, a general structure for the 
A-type superpotential was proposed, which had the integrality properties that 
we referred to above. 

It is of interest to ask the question whether the A-brane world-volume su- 
perpotential, together with the open-string instanton corrections, can be directly 
computed without explicit recourse to mirror symmetry, particularly in the frame- 
work of the GLSM |TB[ . In this paper, we will take the first steps in this direction. 
We will show that the description of the A-branes in the linear sigma model frame- 
work suggests a natural and concrete description of these (partially compactified) 
open-string instanton moduli spaces. These instanton moduli spaces will be seen 
to be intimately related to some auxiliary toric varieties that we shall refer to as 
boundary toric varieties. For these boundary toric varieties we will write down 
the appropriate Gelfand-Kapranov-Zelevinsky (GKZ) differential equations, by 
an extension of the methods developed for the closed string case, particularly 
for non-compact CY manifolds by [|i~7| ,|l8|l. We will then show how the relevant 
solutions of the GKZ equations describe the world- volume superpotential for the 
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A-branes in the examples considered by |||9[]. We will also show how to obtain 
the analogue of the open-string mirror map in this context, using in part, an 
ansatz due to Aganagic, Klemm and Vafa. 

This paper is organised as follows. In section 2, we describe some basic facts 
about the Gauged Linear Sigma Model. After describing the construction of A 
and B-type superspace, we discuss the boundary conditions of || that describes 
A-branes. In section 3, we discuss the GLSM for A-branes and the boundary 
conditions therein. In section 4, we discuss the topological aspects of the A- 
model, before moving on, in section 5, to the open string instanton moduli space. 
After discussing briefly certain issues of stability, we introduce the boundary toric 
variety. In section 6, we consider aspects of local mirror symmetry and the GKZ 
system of equations for our boundary toric variety. We discuss how to obtain the 
bulk and boundary periods that are relevant for the computation of the A-brane 
superpotential in the mirror. Section 7 deals with some examples in order to 
illustrate the methods. We conclude with some observations in section 8. 



2 Background 

The supersymmetry generators for (2, 2) supersymmetry can be written as (in 
the notation of [16| ) 

&susy £ Qa Qa£ 

= -e_Q+ + e+Q_ + e_Q+ - e+Q_ (2.1) 

In superspace with coordinates (x m ,9 ± ,9 ± ), the supersymmetry generators have 
the following representation 

Q± = ^+i6 ± {d ±d 1 ) 

Q± = -~-i9 ± (d ±d 1 ) (2.2) 

On the boundary, one set of linear combinations of this supersymmetry is 
preserved. Let us denote the unbroken supersymmetry generators by Q and Q. 
Let us denote the boundary in superspace by the coordinates (x°,9, 6). Then, 
the supersymmetry generators can be given the following representation 

Q = de +m ° ' ® = ~W~ m ° ' (2 ' 3) 
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We will now work out the precise relationship between the boundary coordinates/ 
super symmetry generators. We will deal with the A- and B-type cases separately. 



A-type superspace 

Under A-type super symmetry, let 

e A = -V2e_ = -V2r]e + , 

where we use the subscript A to denote that we are considering boundaries that 
preserve A-type supersymmetry. Then, on using 

(5 a = ^aQa — ^aQa 

= -e + (Q_ + V Q + )+e + (Q^+ V Q + ) , (2.4) 

we obtain 

Q A = (9±±^Z^ and e A = y/26 + = -V2 V e- • 
B-type superspace 

In a similar fashion one chooses e# = V2e- = \f2r]e + . 

x f-Q+ + yQ-\ _ f-Q+ + vQ- 

Ob = £B 7EF - ZB 



V2 J V V2 
= £bQb ~ £bQb , (2.5) 

from which we identify 

Q B = ~ Q+ V 1 ®" and Q B = -V29 + = v^fl - . 
v2 

2.1 Boundary multiplets 

In boundary superspace with coordinates (x°,9,9), the superderivatives are 

D = oe~ m ° ' D = ~W + m ° • (2,6) 

We will deal with two kinds of boundary multiplets: a boundary chiral multiplet 
and an unconstrained complex multiplet. 

A boundary chiral multiplet satisfying the constraint 

DL = , 
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which has the following component expansion 

L = l + i6\- ieed i . 

An unconstrained complex superfield U with component expansion 

U = u + Gv - 9£ + 99 (D' + iR) (2.7) 

where D' and R are real fields. Under supersymmetry, the transformations are 

5u = — e£ + ev 

5v = -e(-id u + D' + iR) (2.8) 
5£ = -e{id u + D' + iR) 
S(D' + iR) = — ied v — ied ^ 

Now consider a complex superfield with the following gauge invariance 

U^U + L (2.9) 

where L is a chiral superfield. A gauge-invariant observable is given by S = DU . 
The invariance trivially follows from the chirality of L. S has the following 
superfield expansion 

E = DU = £ + 9(D' + iR + id u) - i66d t , (2.10) 

We will identify the real part of U with the bosonic u discussed earlier. Keep- 
ing this in mind, we can partially gauge-fix the imaginary part of U as follows 

lm(U) = 66(R) . (2.11) 

Note that the gauge choice is algebraic and hence we will call it the Wess-Zumino 
(WZ) gauge. This leaves us with the following gauge invariance which acts as 

5Re(u) = Re(Z) 

5R=-d Re(l) (2.12) 
In the WZ gauge, the field U has the following superfield expansion 

u = u + ez-ez + ee(D' + iR) (2.13) 

where u is now a real field. Thus, if R = 0, it would be a real superfield. 
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The supersymmetry and gauge transformations which preserve the WZ gauge 

are 

5u = -e£ + e£ + Z 

5£ = -e(id u + D' + iR) (2.14) 
5D' = —ied t; — ied l; 
5R = -d l 

where / is now a real parameter. Notice that R is invariant under supersymmetry 
and transforms as — d l under gauge transformations and thus behaves as a gauge 
field. 

2.2 Decomposition of bulk multiplets: A-type 

First, consider a chiral superfield $ which has the following component expansion: 

$ = $ + y/26 + i/j + + V29-^_ + 26 + 9-F -i9 + 9 + (D + D^cf) 

- i9-6~-(D - D{)<t> + i9~9 + V2Qa(j) - %9 + 9~\f 7 iQa^ 

- iV29 + 9~9~- [(D - Dx)ip+ + V2Qaip-} 

+ iV29 + 9-9 + l(D + D 1 )^ + V2Qaij + } (2.15) 

In order to obtain the A-type superfield, we need to set 9 = \^29 + = —\[2r\9~ 
Thus, we obtain 

$' = $\e+ = -r,e- = + H+ ~ VH- ~ 00(r]F + iDrf) (2.16) 

where the prime indicates an A-type superfield. This is a complex unconstrained 
superfield. 

Next, let us consider a twisted chiral superfield X which has the following 
component expansion 

S = a + iV29 + \ + -iV29~-\- + V29 + 9-(D-iv i) 

- i9~9-{d Q - di)a - i9 + 9 + (d + dja + V29~9 + 9-(d - di)A+ 

+ V29 + 9~9 + (d + 9i)A_ - 9 + 9-9~9 + (d 2 - d^a. (2.17) 

We define the A-type superfield £' as in the case of the chiral superfield. We 
obtain 

£' = S| e+= _^-- = a + ifl(A+ + ??A_) - i99d a (2.18) 
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which is the superfield expansion for an A-type chiral superfield. In addition, one 
obtains a Fermi chiral superfield A whose lowest component is A = (A + — i]X_) 
satisfying 

DA = 2diE' , (2.19) 

where one treats d{E' as an A-type chiral superfield with lowest component d\a. 
The component expansion for A is 

A = (A + -r/A_)+zv / 2#(L>-woi)-2£<9i(T+# [-id (X + - r?A_) + 2idi(A+ + i]X-)] 

(2.20) 

Something interesting happens if D — ivoi is set to zero on the boundary (part 
of which we know is required by the vanishing of ordinary variations). Then, 
using the bulk equations of motion for the A's, one can see that A becomes an 
anti-chiral boundary superfield. 



2.3 The AV boundary conditions 

Aganagic and Vafa || have constructed a class of special Lagrangian(sL) cycles 
in toric varieties by adapting the construction of Harvey and Lawson in C n . We 
shall describe this construction in the language of the GLSM. Let $j (i — 1, . . . , n) 
be the chiral superfields with U(l) charges Q a { which appear in the gauged linear 
sigma model associated with the toric variety. The toric variety is obtained by 
imposing the (bulk) D-term constraint given by 

f^=r a -^Q^| 2 = (2.21) 

i 

where a = 1, . . . , h 1A . After taking into account the gauge invariances involved, 
one ends up with some weighted projective space (actually, a line bundle on that 
space) of dimension {n — h\^). This is the standard procedure called symplectic 
reduction. The Calabi-Yau condition is given by ^ Q1 = 0. 

Lagrangian submanifolds in C n are specified by the following D-term like 
conditions 

Y,q?\&? = c a a = l,...,r (2.22) 

i 

Note however that there is no associated local U(l) invariance as occurs in sym- 
plectic reduction. One further imposes (n — r) additional conditions on linear 
combinations of the phases Im log 0j of the fields 4>i 

^wflmlog0i = O p = l,...,(n-r) (2.23) 
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The (q a , v") specify a Lagrangian submanifold of C n (with the standard symplec- 
tic structure) if 

£of = o . 

i 

In order for this Lagrangian submanifold in C n to descend to the toric variety, we 
choose hi t i of the conditions given above to be the D-term conditions associated 
with the toric variety. The conditions ^ Qt^i — imply that the conditions on 
the phases are invariant under £7(1) gauge transformations. Further, it follows 
that the submanifold is special Lagrangian if 

£tf = ° • 

i 

This is clearly compatible with the Calabi-Yau condition on the bulk £7(1) charges. 

For generic values of c a , the sL submanifolds have a boundary. One can ob- 
tain sL submanifolds without boundary by a doubling procedure. However, there 
exist non-generic situations (such as the case when one of the c a vanish) where 
the doubling is not necessary. In other words, when we take an sL submanifold 
without boundary for generic values of c a to the non-generic value, the sL mani- 
fold splits into twcQ. We shall refer to one component as an half-brane following 
AV ||. In the toric description, half-branes lie on the edges of toric skeletons. 

In order to be more concrete, let us consider a specific example furnished 
by the non-compact Calabi-Yau manifold given by the total space of the line- 
bundle 0(— 3) over P 2 . This non-compact Calabi-Yau manifold is described in 
the GLSM by four fields with £7(1) charge Q = (—3, 1, 1, 1) and the boundary 
conditions of the special Lagrangian submanifold are given by q 1 = (—1, 1,0,0), 
q 2 = (—1, 0, 1, 0) and v = (1, 1, 1, 1). For generic values of c a , the topology of the 
sL submanifolds (after doubling) is M x S 1 x S 1 . However, for suitable choices of 
c a (see below), one obtains half-branes which are topologically C x S 1 , 

Following the discussion in AKV ||, there are three phases (half-branes) of 
interest: Phase I given by c 2 = and < c 1 < r, Phase II given by c 1 = and 
< c 2 < r and Phase III is given by c 1 = c 2 and c 1 < 0. The relevant phase 
diagram is as follows 

In phase I, for large values of c 1 and r when quantum corrections are small, c 1 
measures the size of a holomorphic disc ending on the half-brane. This is the real 
part of a complex modulus, w 1 , whose imaginary part corresponds to turning 
5 or a suitable power of two depending on the number of doublings involved 



s 




Figure 1: Phases (half-branes) for 0(— 3) on P 2 



on a Wilson line wrapping the S 1 of the half-brane. w 2 is the other modulus 
associated with the other S 1 which appears for generic values of c a vanishes on 
the half-brane. These values suffer quantum corrections which become relevant 
for small values of c 1 and r. We will denote the quantum corrected quantities by 
unhatted quantities. However, w 2 will not be independent and will be determined 
as a function of w 1 on the half-brane. 



2.4 The AV/AKV computation of the A-brane Superpo- 
tential 

As we mentioned in the introduction, the essential idea of [§,[| is to compute the 
superpotential for the A-branes described in the last subsection by using mirror 
symmetry to calculate this in the B-model. Here, we will briefly review their 
arguments. We will be interested in the case of a space-filling D6-brane wrapping 
a sL cycle of a non-compact CY three-fold. 

It was shown in [12] that the amplitude of interest, the (instanton generated) 



disc amplitude in the A-model, is of the form 



°° d - - 
W = J2J2 ^q nk s nrh (2.24) 



71=0 m,k 



where s a = e " ,Q , w a is the (complexified) volume of the holomorphic disc with 
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boundary ending on a particular one-cycle on the sL submanifold and q a = e~ ta , 
t a being the closed string Kahler class, m are the winding numbers associated 
with the holomorphic disc i.e., an element of H\ of the sL submanifold and k is 
an element of H2 of the CY. The integer n relates to the multiple coverings of 
these. Further, this prediction involves the highly non-trivial integers £, that 
count the number of primitive disc instantons in this homology class. These have 
the interpretation of counting the degeneracy of domain wall D4-branes ending 



on the space-filling D6-brane |12| , which wrap a 2-cycle in the CY (whose Kahler 
class is given by k) and m denotes the wrapping number around the boundary. 
Since a direct computation of the superpotential in the A-model involves 



summing up open-string instantons, AV and AKV use the mirror map of [1C] to 
map this computation to one involving holomorphic Chern-Simons theory in the 
mirror B-model. For the case when the B-brane is of complex dimension one, 
and wraps a Riemann surface C in the mirror CY, the dimensional reduction of 
the Chern-Simons theory on C was computed in || in terms of the holomorphic 
three-form of the Chern-Simons theory restricted to the brane and two scalar 
fields (denoted by u and v parametrising normal deformations of the B-brane 
in the mirror CYQ The final result for the superpotential is given by the the 
Abel-Jacobi map associated with the one-form vdu 

nU 

W(u) = / v{u)du (2.25) 

J u* 

where u* is the 'location' of the half-brane. v is fixed on the B-model side by the 
condition that (u, v) be points on the Riemann surface C. However, from the 
A-model viewpoint, it suffices to note that 

v = d u W{u) (2.26) 

which is obtained from the Abel-Jacobi map. 

There are a couple of issues that one has to address before one can translate 
these results in the B-model to the A-model side. First, there is an SL(2, Z) 
ambiguity in the parametrisation of the B-model variables given by 



(2.27) 



h\ fa b\ ( u 
y]_ \c d) \v 
respectively, in the notation of the previous subsection. For the 
convenience of the reader, we shall use the convention of AKV in this subsection. 
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Thus, the superpotential given above is also dependent on the choice of u and v. 
Second, the variables which appear in the B-model side are not the same as the 
ones that appear in the A-model. The change of variables that relates the two is 
the open-string mirror map. 

However, in the limit where quantum corrections are suppressed in the A- 
model, i.e., where the volumes of the holomorphic disc and the CY are large, 
one can relate the two sets of variables. However, the condition that Re(tt) be 
the volume of the holomorphic disc and v — > in this limit, breaks the SL(2, Z) 
down to the subgroup given by 



u\ / 1 n\ ( u 
WVO 1 J \v 



(2.28) 



Thus, imposing the classical limit leaves us with an ambiguity parametrised by 
the integer n. It turns out that for cases where the A-model is dual to large N 



Chern-Simons theory |L2], [H], [TJJ , this integer reflects a framing ambiguity in the 
Chern-Simons theory. 

The framing ambiguity has recently been studied in the context of the Wilson 
loop variables of Chern-Simons theory in |S| , and evidence has been found that 
the above mentioned duality is more appropriate in the context of U(N) Chern- 
Simons theories rather than the SU(N) ones originally considered in the proposal 



of [0. Initial works on Knot invariants were insensitive to this because the 
computations were performed in the standard framing, where the U(N) and 
SU(N) Knot invariants are the same. However, hints of this had already appeared 
in ||24)| . In explicit computations of certain U(N) Wilson loop operators have 
been carried out, and the results have been shown to conform to the integrality 
predictions of RI5| . 

Aside from the ambiguity mentioned above, one has to obtain the open-string 
mirror map in order to obtain the A-model results. It has been argued in AKV 
that the appropriate variables are given by || 

u = u + A(q) (2.29) 
dW 

* - -ar (2-30) 

where u receives corrections only from closed-string instantons. u is computed by 
identifying it with the tension of the domain wall which is mirror to the D4-brane 
in the A-model. The variables u and v are conjugate fields of the holomorphic 
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Chern-Simons theory. The second equation corresponds to the condition that v 
remains conjugate to u in the quantum theory. 

As we will see in the next section, the B-model variables, i.e., the unhatted 
variables, naturally appear in the GLSM that we construct to describe these 
half-branes. 

3 Linear Sigma Models for A-branes 
3.1 The basic boundary condition 

The basic A-type boundary conditions that one considers in a GLSM with n 
chiral superfields and r U(l) vector multiplets are 

5 = ^-E^| 2 = a = !,•••»'■ (3-1) 

i 

The bosonic boundary conditions (in the matter sector) are completed by the 
following (n — r) conditions 

^^Imlog0 4 = O b = l,...,(n-r) (3.2) 

i 

where the Lagrangian condition implies that the u's satisfy 

i 

This condition permits us to identify, in the toric description, the v\ with the 
integral generators of the one- dimensional cones in the fan A, the combinatoric 
data associated with the toric variety. 

The boundary conditions on the vector multiplet are given by the superfield 
condition 

K = (3-3) 

The easy way to understand this is to notice that in the topological A-model, the 
field a is proportional to the Kahler form u whose restriction on a Lagrangian 
submanifold is zero. The vanishing of the boundary terms in the ordinary varia- 
tion is obtained by the conditions 

«oi = (3.4) 
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The above set of boundary conditions form a consistent set of boundary conditions 
in the GLSM which also have a proper NLSM limit in the sense used in [ 19 , pOf . 
Let us denote by \B)o the boundary state (in the GLSM) that corresponds to 
these boundary conditions. 

This reference state plays a role similar to the D6-brane in the context of 
B-branes in the GLSM pl|,|22|. It is not hard to see that it corresponds to the 
D6-brane in the mirror Calabi-Yau manifold. It is also rigid i.e., there are no 
moduli associated with this sL submanifold. 



3.2 The AV boundary conditions in the GLSM 



The boundary conditions that AV choose can be implemented in the GLSM in a 
simple fashion^. In the matter sector, it corresponds to simply replacing some of 
the v 's by boundary conditions of the form| 
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In the GLSM, we will realise these boundary conditions as low-energy conditions 
using complex boundary supermultiplets U a . The corresponding boundary state 
can be represented by the path-integral 



AV 



[DU] exp^ + S 2 + S 3 ) \B) { 



where 



u a + u a 



-^^Imlog^ 



This action is invariant under the gauge-invariances 



(3.5) 



(3.6) 



Im log $^ 



U a + L a 
Imlog$. + qf 



L a + L c 



(3.7) 
(3.8) 



where L Q are chiral superfields. This gauge invariances enable us to identify the 
real parts of the lowest components of the superfields U a and with "coordinates" 



7 A similar construction was pursued in Ea| though it differs from ours in some details. 
8 The angles dual these are given by g^Imlog^i, where are the dual basis vectors (in C n ) 
to the qf i.e., they satisfy: ~q\qf = 6%, £ 4 qj^Q" = and g>f = 0. 



13 



on an appropriate three-cycle. The action Si is also invariant under bulk gauge 
transformations. On using the fact that D a are generators of the bulk gauge 
invariance, one obtains that 

D a \B) AV = , 

given that S2 and S3 (given below) are also gauge- invariant. 
The kinetic energy for E a = DU a 

S 2 = ~J dx°d 2 9 E a E a 

= ^J dx ° (%«doZa + \{D' a f + ^(R a + d u a ) 2 ^j (3.9) 

where we have introduced the coupling constant e& on the boundary. This has 
mass dimension half. 

We can now introduce the boundary analogue of the F.I. term: 

S 3 = -Re (w a j dx°d9 E a ^j = - J dx° (c a D' a - ^-(i? a + d u a ^J (3.10) 
where w a = c a + i^- is a complex parameter. 



3.3 Taking — > 00 

In the limit — > 00, the fields in the multiplet E a impose the constraints 

'52q?@ i 'iI> + i + vM-i)=0 (3.11) 

i 
i 

V Yl tf& F i + + 1 E ^(0^10, - = • 

i i 

The last three equations are the supersymmetric variation of the first equation. 

The 6 a term is a topological term as one expects. Thus, we can see that the 
real part of R behaves like the bulk D term while the imaginary part should 
behave like v 01 . In order for this correspondence to work, we should require that 
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(after Wick rotation to Euclidean space and x° becomes a circle variable which 
we will call x 2 with periodicity 2tt) 



/ — (R a + d 2 u a ) = -m a 
J 

where m a are some integers. 

4 The topological version of the GLSM 

In order to discuss the topological aspects of the A-model, we need to Wick 
rotate: x° — > — ix 2 (Hence one has d — > id 2 , R — > iR and i>oi "~ > ~ ^12) • Let us 
choose to twist the model so that Q is the BRST operator. For the closed string 
case, the localisation equations are |16| (a = 1, • • • , ^1,1). 




{D 1 + iD 2 )<f) = 
D a + ^i2 = 
Fi = 



(4.1) 



a a = 



Let n a be the instanton numbers defined by 




(4.2) 



The action is then equal to 



2mt a n, 



a 



(4.3) 












(4.4) 



in addition to the boundary condition 




(4.5) 



The winding numbers are defined by 




(4.6) 
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The contribution of the boundary terms to the action is 



SWy = -2nis a m c 



(4.7) 



Thus, the total action in the topological theory is given by 
S = S huik + S hdzy = -2nit a n a - 2nis a m a 



4.1 Observables in the topological theory 

The gauge-invariant BRST closed observables in the closed string sector are a a . 
These have ghost number two. In the open-string sector, the observables in the 
NLSM limit should be in one-to-one correspondence with closed forms on the 
Lagrangian submanifold. In the AV construction, we assume that a subset of 
closed one-forms is given by du a , where we consider u a to be (compact) coordi- 
nates (angles) on the Lagrangian submanifold. Thus, in the topological model, 
the corresponding observables are £ a (which are naively BRST exact). Thus one 
has 





Observable 


Ghost No. 


Bulk 




2 


Boundary 




1 



5 The open-string instanton moduli space 

The open-string instanton moduli space can be defined to be the moduli space 
of the localisation equations that we derived in the previous section. We assume 
that the worldsheet is a disc. The relevant equations are: 



(D 1 +iD 2 )(P l = 
D a + v a l2 = 
D' a + (R a + d 2 u a )\ 9D =0 
J^uflmlog&|a£> = 



(5.1) 
(5.2) 
(5.3) 
(5.4) 



The first two equations are bulk equations while the last two equations are on 
the boundary of the disc. Thus, we are interested in the moduli space of these 
equations modulo bulk as well as boundary (real) gauge invariances. Further, 
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the instanton moduli spaces are fixed by the instanton numbers n a as well as the 
winding numbers m a . 

The first equation is invariant under complex gauge transformations. As in 
the bulk case, we shall consider the last two equations as gauge-fixing of the 
imaginary part of the gauge transformation. Thus, we shall be interested in the 
moduli space of the equations 

(D 1 + iD 2 )(f)i = (5.5) 

modulo complex gauge invariances. Under suitable circumstances i.e., when one 
is dealing with a stable situation, this moduli space is identical to the moduli 



2vm a + e 2 / Q, a |0i| 2 = eVvol(JJ) (5.6) 



space of the full set of equations including (|5.2| ) and ( |5lf ) modulo real gauge 
invariance. In a suitable gauge where t>f = 0, this equation states that (pi is a 
holomorphic function on the disc. In the closed string case, this implies that (pi 
are sections of 0(rii = J2 a Q°:n a ) in the sector with instanton number n a . 



5.1 Stability 

In order to obtain the conditions under which stability is obtained (see for instance 
pj|j), we shall first consider the bulk equation (|5.2|). Let us choose the worldsheet 
be a disc D. Integrating (|5.2|) over the disc we obtain 

,2 / n»l A |2 2_a 

JD 

Let us consider the case when e 2 r a ^> (for some fixed a) and all fields (pi whose 
£7(1) charge is negative are set to zero. One then obtains the following inequality 

, \ 2?rn a . . 

™m > — (5.7) 

It is not hard to see that for fixed n a > 0, this condition can always be satisfied for 
large enough e 2 r a . Thus stability is guaranteed in this case. A similar argument 
occurs for the boundary equation ( p.3| ) leads to the condition 

vol(dD) > (5.8) 

where stability is again guaranteed for large enough e\c a for m a > 0. The 
conditions involving the phases (eqn. ( |5.4| )) can always be imposed and hence do 
not play a role in considerations of stability. 
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Let us now consider the situation where the disc is large i.e., can be ap- 
proximated by the complex plane. Then, one has the following interpretation 
for the bulk instanton numbers and the boundary winding numbers. Consider 
rii = Q1n a + q"m a . Then, one has (as solutions of the localisation equations) 

1. (pi = whenever rij < 0. 

2- (pi = Y^jLo&hj ^ otherwise 
where z is a (complex) coordinate on the disc. The winding numbers (instanton 
numbers) fix the leading power of z (number of zeros on the disc). Let us denote 
this space by Y n a^ ma . 

Thus, the moduli space is given by the space of (pij for all rij > modulo 
gauge invariances. This is fixed by 

1. The bulk gauge invariance is fixed by a D-term of the form 

i \j=0 / 

and the associated U(l) gauge invariance. 

2. The boundary condition on the disc, eqn. ( |5.3|) , becomes a constraint of 
the form 

i \j=0 / 

3. The boundary conditions in ( |5.4j ) can be considered to arise from the gauge 
fixing of the £7(1) gauge invariance 

4>i,j — > e lv ' 9 4>i,j ■ 

In this case, the moduli space is the variety given by 

■M-n a ,m a = (Yn a ,m a ~ Fn a ,ra a ) jT (5-9) 

where T acts as 

fa -> A Q ? fa for A G C* , (5.10) 

<j>i,j -> \ q? (f>ij for A G R + (5.11) 

-> e iv * e <p id for real 9 . (5.12) 

F na ,m a is the set of points in Y natTria for which the D-term constraints given above 
are not satisfied. Note that in the zero instanton sector (i.e., n a = 0,m Q = 0), 
one obtains the sL submanifold as expected. 
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5.2 The boundary toric variety 

Suppose, we treat the D-term-like boundary conditions which are specified by gf 
on par with the bulk D-terms specified by Q". The special Lagrangian condition 
Yli Qi = now becomes the Calabi-Yau condition on a toric- variety which we 
shall call the boundary toric variety. The boundary moduli now become bulk 
moduli on the toric variety. One can see that the spaces Y na>ma and M. na ,m a 
obtained in the previous subsection play a similar role for the boundary toric 
variety. A well known method of computing mirror periods of a Calabi-Yau 
hypersurface is by the Picard-Fuchs system of differential equations. By treating 
the bulk and the boundary D-term constraints on the same footing, one might 
hope to compute the period integrals of the boundary toric variety in a similar 
fashion. This should, in principle, directly give the necessary variables for the 
computation of the worldvolume superpotential for A-branes. In the following 
sections, we will proceed with such a computation, and show that the boundary 
toric variety indeed reproduces the structure of the A-model amplitude predicted 
in 0. 

6 Summing up open-string instantons 

6.1 The GKZ equation for the bulk GLSM 

As discussed earlier, the toric data is seen in the GLSM construction through the 
U{\) charges, (a = 1, . . . , r), of the chiral superfields and v 1 ? (b = 1, . . . , (n—r)) 
which are the integral generators of one-dimensional cones in the fan A of the 
toric variety. Consider the differential operators constructed from this data 

^ - n (if -n(£f ? («) 

Q a >0 \oatS Qt<Q \oa % j 
n f) 

i=0 1 

with (3 = (—1, 0, . . . , 0) and are coordinates on C n . The two operators given 
above define a consistent system of differential equations, due to Gelfand, Kapra- 
nov and Zelevinsky, the GKZ hypergeometric system with exponent (3 

Cqo. 7r(a) = , Z b 7f(a) = (6.3) 
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It is known that the periods (associated with three-cycles) of the mirror 
Calabi-Yau manifold are solutions of a Picard-Fuchs differential equation which 
naturally arises in the special geometry associated with vector multiplets with 
N = 2 supersymmetry in four dimensions. It turns out that for the CY mani- 
folds which arise as hypersurfaces in toric varieties, the GKZ equations described 
above are closely related to these Picard-Fuchs equations. In the cases that we 
will be interested in, the above system of equations provide the Picard-Fuchs dif- 
ferential equations. In general however, one may need to supplement the above 
system by additional differential operators to convert the GKZ into the relevant 
Picard-Fuchs equation. 

The second set of differential equations, i.e., those given by the Z b , are solved 
for by a change of variables dj — > z a given by 



and the change 7f(a) = ^Tr(z). Further, it is useful to rewrite the differential 
operators Cq* in terms of logarithmic derivatives O a = z a -J^-. 

Let us assume that z a — > is the large- volume limit. There are two classes 
of solutions (periods) of the GKZ equation which are of interest (see W7\ for the 



case of compact CY manifolds and [18| for the non-compact case [local mirror 
symmetry]). 

1. Those with logarithmic behaviour as z — ► 0. These are important in pro- 
viding the change of variables which form the mirror map. This also relates 
the the GLSM variables with the NLSM (large-volume) variables. This is 
clearly seen in the work of Morrison and Plesser |27J . 

2. Those that behave as (logz) 2 as z — > 0. Using special geometry, these are 
related to the first derivatives of the prepotential. This prepotential encodes 
the corrections due to closed-string instantons to the classical one obtained 
from the triple intersection numbers of two-cycles on the CY manifold. 



6.2 The GKZ system for the LSM with boundary 

We will now show that we can sensibly fit the case of LSM with boundary into this 
framework and write down the corresponding GKZ equations. We will explicitly 
show this for all the cases that have been considered by AV and AKV. Loosely 
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speaking, the AV boundary conditions correspond to replacing some of the i>'s 
with g's. This suggests that the GKZ obtained by the set of charge vectors 
(Qti 1i) might be relevant for obtaining the open-string instanton effects just as 
the bulk GKZ for the closed-string case. Consider the following operators 

- n(£f-n(£f Q? (6 4) 

Qf>0 v l/ Qf<0 v %/ 

^ - n(£f-n(£f («■«) 

i=0 1 

This is the GKZ which is naturally associated with the boundary toric variety. 
As before, we will solve the differential equations given by Zp by a change of 
variables 

i i 

where s a are related to boundary moduli just as the z a are related to bulk moduli. 
It is useful to separate the remaining differential operators into those that arise 
from the Q\ (we will call these £buik) an d those that arise from the qf (we will call 
these £bdry)- Note that £tmik is not quite the differential operator which appears 
in the closed-string GKZ - it also involves the boundary variables. 



6.3 Identifying the solutions of interest 

We will restrict our attention to the cases of the half-branes with topology CxS 1 
and let s = exp(-w) be the modulus associated with the S 1 . The limit where 
quantum corrections are small is given by z, s — > 0. 

By construction, the bulk periods are solutions of the differential equations 
involving £buik- In particular, the solutions that provide the closed- string mirror 
map are solutions with asymptotic behaviour logz as z, s — > 0. The AKV ansatz 
plays a similar role in obtaining the open-string mirror map i.e., we look for 
solutions of £buik with the property 

w = log s + A(z) 

It is not hard to see that the above ansatz is not compatible with £bd ry since A 
is not a function of s. 
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The solution associated with the superpotential is obtained by an ansatz of 
the form 

d Q W(w) = J2 a $,m) z l s m 

k,m 

where we require that the zero instanton contribution vanish i.e., a(0, 0) = and 
the only non-vanishing a(k, m) are as given by the considerations of section 5. To 
be more precise, we set a(k, m) = whenever the corresponding instanton moduli 
space does not exist. Unlike the solution associated with the open-string mirror 
map we require that this solves the full boundary GKZ differential equation. 



7 Examples 

7.1 0(-3) on P 2 
7.1.1 Bulk Periods 

The GKZ equation is given by 

[03 + ^(3e 2 + 2)(30 2 + l)30 2 ] tt = (7.1) 

The basic solutions which we consider are fixed by the behaviour near z = 0: (i) 
the constant solution; (ii) the log ,2 solution and (iii) the (logz) 2 solution. Using 
the Frobenius method, one obtains 

7T = 1 (7.2) 

(nlr n 

n=l v ' 

* = U + C2 + E(") n 7^ (--+ it-)- (7-4) 

where C = Er=iH n (^f- 

The mirror map relates the large- volume coordinate t to 7Ti. Thus, we obtain 

t = -W-o = -log,-B-)" M 3- (7-5) 

n=l ^ ' 

which has the right behaviour z = e~* as z — > and the terms in the summation 
are instanton corrections. The inverse relation is of the form 

z = q + 6q 2 + 9q 3 + 56g 4 - 300g 5 + 3942g 6 - 48412g 7 + • • • (7.6) 

where q = e~ l 
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7.1.2 Phase I 

The GKZ equation for the boundary toric variety is given by 

[9 2 (6 2 + Q s ) + z(3Q z + Q S + 2)(36 2 + Q s + 1)(36 2 + Q a )] ix = (7.7) 

[(6 2 + s )+ S (30 2 + s )] tt = (7.8) 

The first equation represents £tmik and the second represents £bdr y for this exam- 
ple. Since c 2 = in this phase, we have "switched-off" the coordinate associated 
with c 2 . 

The boundary flat coordinate u satisfies the first of the two equations given 
above. For z, s — * 0, one requires that u = — logs. It can be argued that the 
general solution should take the form 

u = — logs + A(z) 

i.e., the boundary flat coordinate differs from the GLSM coordinate only by closed 
moduli. One can then easily prove that 

u — — log s — - (t — i) + in 

is a solution of the first of the two equations (i.e., £buik) that make up the 
boundary GKZ equations. 

Defining s = e~ u , one obtains the following 

s = -5(1 + 2q - q 2 + 20q 3 - 177g 4 + 1980g 5 - 24023g 6 + 97684g 7 + • • •) (7.9) 

This equation along with the earlier equation for z, gives us the required change 
of variables from z, s to the flat coordinates s and q. 

We are now in a position to compute the boundary superpotential. We require 
that it be a solution v of both equations of the boundary GKZ. It should satisfy 
the following conditions: (a) v should vanish and be regular at z = s = 0. (b) 
By consider the instantons in the GLSM, one obtains the vector kQ + mq l = 
(—3k — m, k + m,k, k). In phase (i), we need both 0! and 3 to be non-vanishing. 
This is achieved if k > and k + m > 0. This suggests the following ansatz for 
this solution: 

oo oo 
k=0 m=—k 
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m\k 





1 


2 


3 


4 


5 


-5 


o 


o 


o 


o 


o 


5 


-4 














-2 


28 


-3 











1 


-10 


102 


-2 


o 


o 


-1 


4 


-32 


326 


-1 





1 


-2 


12 


-104 


1085 


1 


1 


-1 


5 


-40 


399 


-4524 


2 





-1 


7 


-61 


648 


-7661 


3 





-1 


9 


-93 


1070 


-13257 


4 





-1 


12 


-140 


1750 


-22955 


5 





-1 


15 


-206 


2821 


-39315 


6 





-1 


19 


-296 


4450 


-66213 


7 





-1 


23 


-416 


6868 


-109367 



Table 1: Degeneracies for phase I of 0{— 3) on P 2 



It may seem that this solution cannot be regular at z = s = 0. But, since we 
require that c 1 < r, this requires the limit be take such that z/s — > 0. Thus z k s m 
vanishes for z = s = as long as k + m > 0. 
The solution is given by 

v = V V (-l) k+m r (3fc + m) z k s m (7 10) 

Putting in the expressions for the flat coordinates, by making the substitutions 
as in (|7.6|) and ( |7.9|) , and comparing with the expression 

v = -J2^d k , m \og(l-q k s m ) = ^d k , m (q k s m ) n . (7.11) 

k,m k,m,n 

Table 1 lists the degeneracies (4 jm of disc instantons. 



7.1.3 Phase III 

The charges relevant for this phase are: Q = (—3, 1, 1, 1) and q = —q 2 = 
(1, 0, — 1, 0). The relevant GKZ equations are 

[e 2 (e 2 - e r ) + z (3e 2 - e r + 2) (3e 2 - e P + 1) (3e 2 - e r )] n = 

[(36 2 - r ) + r (8, - 8 r )] vr = (7.12) 
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where the boundary modulus is represented by r in this phase. In this case, as 
before, the boundary flat coordinate v has a general solution 

v = log r+ \(t-i) + m (7.13) 
3 

and defining f = e v , we have the expansion 

r = -f [1 - 2q + 5q 2 - 32q 3 + 286g 4 - 3038g 5 + 35870g 6 + • ■ ■] (7.14) 

The allowed instanton numbers in the GLSM are given by looking at kQ + mq = 
(m — 3k, k, k — m, k). Since <fto ^ in this phase, we need m — 3k > in addition 
to k > which suggests an ansatz of the form (for the superpotential) 



oo oo 
fc=0 m=3k 



' = Yl a ( k > m ) * k r ™ ( 7 - 15 ) 



As can be easily seen, the GKZ equations ( |7.12| ) are solved with the choice of 

a(k, m) as 



a(k, m) = {-) k+m - ^ m k \ 7 (7.16) 

v ' 1 K ' r(l -3k + m)T(k + 1) 2 V 1 

The degeneracies that are obtained from these can be seen to be consistent with 

those obtained by AKV (see Table 2). 

7.2 0{-l) + 0(-l) on P 1 

The boundary toric variety is characterised by the charges Q = (1,1,-1,-1), 
and q = (0, 0, 1, — 1). Let us consider the phase where c 1 = and c 2 > (phase II 
of AV). The relevant GKZ equations can easily be seen to be 

[o 2 z -z(e z + e s )(e z -e s )]it = o 

[(e s -e 2 ) + s (e s + e 2 )]7r = o (7.17) 

Where, as before, z = e~* and we define s — e v In this example, there are no 
corrections to the open string moduli coming from closed string ones, and one 
can simply write 

v = v + in , u = u + in (7-18) 

so that s = — s. In this case, the allowed instanton numbers are given by the 
vector kQ + mq = (k, k, —k + m, —k — m), and since 03 7^ 0, we would require 
m > k, which suggests an ansatz of the form 



00 00 

u = 

k=0 m=k 



J2J2a{k,m)z k s m (7.19) 
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m\k 





1 


2 


3 


4 


5 


3 





-1 


3 


-18 


153 


-1560 


4 





-1 


4 


-20 


160 


-1595 


5 





-1 


5 


-26 


196 


-1875 


6 





-1 


7 


-36 


260 


-2403 


7 





-1 


9 


-52 


365 


-3254 


8 





-1 


12 


-76 


528 


-4578 


9 





-1 


15 


-111 


784 


-6627 


10 





-1 


19 


-161 


1176 


-9800 


11 





-1 


23 


-230 


1777 


-14720 


12 





-1 


28 


-323 


2678 


-22347 


13 





-1 


33 


-446 


4018 


-34141 


14 





-1 


39 


-605 


5968 


-52289 



Table 2: Degeneracies for phase III of 0(— 3) on P 2 



The appropriate a(k,m) are easy to find, and, defining z = e *, the solution of 
the GKZ equations can be seen to be 

, k+m T(k + m) 



u 



T(m- k + l)T(k+ l) 2 



z k s m 



(7.20) 



fc=0 m=k 

We now reproduce some of the degeneracies c4, m for this example in Table 3. 
The numbers agree with those in || up to an overall sign. 



7.3 Degenerate Limit of P 1 x P 1 

In this subsection, we will calculate the superpotential for another example, that 
of a non-compact Calabi-Yau, containing P 1 xP 1 . We will consider the theory with 
one of the closed string moduli goes to infinity, and will consider the phase where 
the relevant charges that defines the other closed string modulus and Lagrangian 
submanifold are Q = (—2, 1, 1, 0, 0), q = (—1, 0, 0, 1, 0). The GKZ equations for 
the boundary toric variety defined by these charges are 

[9 2 - z (26, + 6 S + 1) (29, + e.)] vr = 

[<d s + s (2<d z + 9 S )] vr = (7.21) 
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m\k 





1 


2 


3 


4 


5 


11 





1 


-30 


390 


-2730 


11466 


12 





1 


-36 


556 


-4690 


24024 


13 





1 


-42 


770 


-7700 


47124 


14 





1 


-49 


1040 


-12152 


87516 


15 





1 


-56 


1375 


-18564 


155195 


16 





1 


-64 


1785 


-27552 


264537 


17 





1 


-72 


2280 


-39900 


435708 


18 





1 


-81 


2871 


-56520 


696388 


19 





1 


-90 


3570 


-78540 


1083852 


20 





1 


-100 


4389 


-107255 


1647455 


21 





1 


-110 


5341 


-144210 


2451570 


22 





1 


-121 


6440 


-191180 


3579030 



Table 3: Degeneracies for phase II of 0{-l) + £>(-!) on P 1 



where we have defined z = e 1 and s = e v . The mirror map dictates that the 
boundary flat coordinate in this case is given by 

s = -(l + q)s (7.22) 

where s = e v , and similarly for the other open string modulus, q = e~ f is defined 
through 

q (7.23) 



The allowed instanton numbers are determined from the vector kQ + mq = 
(—2k— m, k, k, m) and since 0i, 2 , 04 are non- vanishing, we require that k, m > 0. 
Writing an ansatz for the solution to the GKZ equations, 



oo 

// 

fc=0,m=0 



= a(k,m)z k s m (7.24) 



The coefficients are easily determined to be 
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m\k 





1 


2 


3 


4 


5 


11 





1 


36 


676 


8281 


74529 


12 





1 


42 


920 


12936 


132496 


13 





1 


49 


1225 


19600 


226576 


14 





1 


56 


1600 


28896 


374544 


15 





1 


64 


2055 


41616 


600935 


16 





1 


72 


2601 


58680 


938961 


17 





1 


81 


3249 


81225 


1432809 


18 





1 


90 


4011 


110550 


2140369 


19 





1 


100 


4900 


148225 


3136441 


20 





1 


110 


5929 


196020 


4516475 


21 





1 


121 


7112 


256036 


6400900 


22 





1 


132 


8464 


330616 


8940100 



Table 4: Degeneracies for the degeneration of P 1 x P 



Substituting in ( fi 7 .24| ), with the summation excluding (0, 0), and, substituting for 
z and s in terms of q and s, we recover the integers c4 iJn which are consistent 
with the results of AKV (see Table 4). 



8 Conclusion 

In this paper, we have constructed a GLSM with boundary which corresponds 
to D-branes associated with the class of special Lagrangian cycles proposed by 
Aganagic and Vafa. The considerations of the topological field theory associated 
with the GLSM naturally leads to the an associated boundary toric variety. We 
have seen that the GKZ associated with the boundary toric variety provides us 
with the open-string mirror map and enables us to verify the integrality conjec- 
tures for the degeneracies associated with open-string instantons. We also obtain 
natural models for the open-string instanton moduli spaces. We have however 
not directly computed the degeneracies in the topological model - this involves 
a generalisation to the GLSM with boundary, of the methods used by Morrison 
and Plesser for the closed-string case |27| . The methods used in this paper clearly 
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have a generalisation to cases involving compact CY manifolds i.e., those that 
include a bulk superpotential in the GLSM. However, there might be cases where 
the moduli are related to non-toric deformations for which special techniques are 
necessary. 

The framing ambiguity is quite similar to what may happen in the case of a 
conifold singularity, in say the quintic. In that case, one has an S 3 which shrinks 
to zero size. Its dual three-cycle is a T 3 . The monodromy action around the 
conifold singularity generates the same subgroup of SL(2, Z) that we saw in the 
framing ambiguity. The volume of the T 3 thus suffers from a similar ambiguity in 
the sense that one can always add an integer times the S 3 volume (which clearly 
vanishes at the conifold singularity) to this volume. However, in the case of the 
quintic, this ambiguity can be fixed by global considerations. In a similar fashion, 
we expect the framing ambiguity to be associated with a monodromy of the GKZ 
solutions and that in the compact case, this ambiguity should be lifted by global 
considerations. 

A 'derivation' of the GKZ equations that we use can be obtained in a fashion 



similar to the closed-string case. By extending the Hori-Vafa |I(J mirror map, one 
can represent the periods on the B-model side as a functional integral over some 
dual fields with appropriate superpotentials corresponding to all the D-terms, 
bulk and boundary, that appear on the A-model side. Using this representation, 
it is possible to derive the appropriate GKZ in a manner similar to that of Hori 
and Vafa for the closed-string periods. From a somewhat different view point, 
one notices that the periods can be represented as an integral of the holomorphic 
three-form over the three-cycle and then one may obtain the differential equation 
satisfied by the periods. The GKZ equations that we use can be derived by 
replacing the integral over the holomorphic three-form by the holomorphic Chern- 
Simons action. 

The AV boundary conditions that we have implemented in the GLSM corre- 
spond to two-branes on the B-model side given by intersections of hyperplanes. 
It is of interest to extend our methods to construct the GLSM for the A-model 
mirror to D-branes associated with vector bundles p2] . 



While this manuscript was being readied for publication, a paper |HJ] appeared 
which in part, overlaps with the contents of section 6 and 7. 
Acknowledgments One of us (S.G.) would like to thank the Theory Division, 
CERN and two of us (S.G. and T.J.) the Abdus Salam ICTP for hospitality while 
this work was being completed. 
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